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The Entrainment-Limited Evolution of FRII Sources: 
Maximum Sizes and A Possible Connection to FRIs 
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ABSTRACT 

We construct a simple theoretical model to investigate how entrainment gradually 
erodes high-speed FRII jets. This process is described by embedding a mixing-layer 
model developed originally to describe FRI objects in a self-similar model for the lobe 
structure of classical FRII sources. Following the classical FRII models, we assume 
that the lobe is dominated by the particles injected from the central jet. The entrain- 
ment produces a boundary shear layer which acts at the interface between the dense 
central jet and the less denser surrounding lobe, and the associated erosion of the jet 
places interesting limits on the maximum size of FRII sources. The model shows that 
this limit depends mainly on the initial bulk velocity of the relativistic jet triggered. 
The bulk velocities of FR lis suggested by our model are in good agreement with that 
obtained from direct pc-scale observations on ordinary radio galaxies and quasars. 
Finally, we discuss how FRIIs may evolve into FRIs upon reaching their maximum, 
entrainment-limited sizes. 
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1 INTRODUCTION 

Radio galaxies appear to come in two fundamentally differ- 
ent types, as en capsulated in the Fanar off and Riley classi- 
fication scheme (Fanaroff fc Rilevlll974l ). FRI sources have 
bright cores and edge-darkened lobes, while FRII sources 
are edge-brightened with hotspots at the end. These mor- 
phological differences suggest that the interactions between 
the radio jets and their environments are very different in 
the two classes, and that their evolutionary tracks may also 
be quite distinct. 

FRII sources are powerful radio sources with fairly 
homogeneous morphologies. They contain highly relativis- 
tic jets extending from the central AGN to very bright 
hotspots surrounded by low surface brightness lobes. Dy- 
namical models for FRII sources are quite successfu l, in- 



dicat i ng that the lobes expand in a self-similar way iFallel 
ll99ll : lKaiser fc Alexandeill 19971 . hereafter KA97). Based on 
this, a range of radio emission models have been devel- 
oped, and these models allow FRII sources t o be tracked 
through the power- linear size (P-D) diagram (iKaiser et al.l 
ll997l : iBlundell et al.|[l999l : iManolakou fc Kirkll20oi v 

FR I sources are much more common than FRIIs 
l|Parma et al.|[200^ ). but they are also more complex and 
have only one common feature: no hotspot at the outer end 
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of the jet. About half of the FRI sources show a fat dou- 
ble morphology similar to FRII lobes, while the rest inflate 
turbulent lobes after passing through a so- called brighten- 
ing point, with plumes or tails at the end (lOwen fc Laind 
1 19891 : lOwen fc Wh"it3 Il99ll : iParma et al.] |2002| ). Modeling 
FRI sources is difficult, as it is ha rd to describe all types 
of FR I sources with a single model. iBicknelj (| 19941 ) tried to 
model the FR I s ources by relativistic conservati on laws, and 
iLaing fc Br"idi3 l|2002ah studied 3C 31 in detail. IWang et al] 
l|2009l, hereafter W0 9 ) ado pted the mixing-layer structure 
from ICanto fc Ragal 1 199ll . hereafter CR91) and built an 
analytical model which could explain the observational be- 
havior of typical tailed FRI sources. 

Generally speaking, FRIIs are more powerful than 
FR Is, with a transition radio luminosity around -Pi7g]y[Hz 
lO^^WHz"^ sr~^. A transition lumino sity also exists in 
the optical band (|Owen fc Ledlowlll994l '). The value of the 
transition luminosity is not precise, as it also depends 
on t he properties of the host galaxies l|Ledlow fc OwenI 
119961) . The origin of the FR I/II dichotomy has been dis- 
cussed extensively in the literature. Studies on compact 
steep-spectrum sources (CSS) suggested that these objects 
are typically young and may general ly evolve into large- 
scale radio sources (iFanti et al.l Il995l ). Among these CSS 
sources, FRIs and FRIIs may have diff erent progenitors 
due to different powers and envir onments l|Alexandedl2000l : 
iKunert-Bairaszewska et al.ll2005l ). However, the transition 
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from FRIIs to FRIs can also possibly take place at a later 
stage of the jet evolution under certain c ircumstance l|Fallel 
ll99ll : lBicknellll995l : lKaiser fc Bestll2007l ). 

The instabilities of jets have been studied by a number 
of numerical simulations, and they find that the jet instabil- 
ity evolution and the large-scale jet m orphology are mainly 
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Meliani et al.|[200^ . Ifowever, these numerical simulations 



more concentrated on the earlier stage of the jet evolutions, 
when the lobe structures around the central relativistic jets 
have not been well developed and differ from that at the 
later jet evolution stage. Meanwhile, the initial set of the 
simulations (e.g. t he radial resolution) may also affect the 
simulation results (|Perucho et al.l [20041 ). Therefore, analyt- 
ical models describing the evolutions of jet stabilities are 
desired for studying the jets at the late stage of their evolu- 
tions. 

Taking the idea of the relativistic mixing layer model 
developed by W09, we can investigate more precisely how 
the central jets of FRIT sources are eroded by the entrain- 
ment due to the interactions with their surrounding lobes. 
More specifically, we describe the entrainment process by 
embedding the W09 mixing- layer model for FRI jets in a 
simple, self-similar model of an FRII radio lobe, and mon- 
itor how the central jet is gradually eroded by the growing 
turbulent shear layer at the interface between the jet and the 
lobe. The goal of this paper is to study if the entrainment 
can play an important role on determining the maximum 
size of an FRII jet at the late stage of its evolution. Our 
basic model for a central jet and shear layer embedded in an 
FR II radio lobe is developed in Section (2] The maximum 
size of FRII jet is calculated and discussed in Section [S] 
We also argue that the resulting dead FR lis will ultimately 
re-emerge as FR I radio sources. A brief sketch of this tran- 
sition process is provided in Section [l] Finally, in Section 
[S] we summarise our conclusions and outline future work 
suggested by our model. 



2 MODEL DEVELOPMENT 

We describe FRII objects by embedding a highly relativis- 
tic central jet inside a surrounding radio-emitting lobe. Al- 
though the lobe density is thought to be very low, we as- 
sume that a turbulent shear layer may nevertheless form at 
the jet-lobe interface. This shear layer will entrain and mix 
material from both regions. This entrainment will gradually 
erode the central jet. More specifically, once all of the highly 
relativistic material in the central jet has been mixed up with 
the lobe material in the shear layer, the central jet is com- 
pletely destroyed. Previous work suggests that the hotspot 
is a very compact, high-pressure region that gives rise to 
strong radio emission (jScheck et al.ll2002l '). The highly rela- 
tivistic central jet may play an important role in energising 
the hotspot as there is large amount of energy injected into 
a small area. Meanwhile, the bulk velocity of the material 
in the shear layer is not as fast as that in the central jet. Al- 
though the shear layer is still supersonic and can form weak 
shocks and working surfaces, it is not powerful or concen- 
trated enough to support a hotspot. Thus, we assume that 
as the central jet is gradually eroded, the hotspot weak- 



Pc Pc 



to 




Lobe 



ti 




Shear Layer 



Refilling flow 







A 

















Figure 1. Sketch of the evolution of a radio outflow. At io, the 
young outflow is showing a FR II morphology. At t\ , the outflow 
is still in FRII phase while the shear layer has already grown. 
At 42, the hotspot vanish and the outflow will transfer into FRI 
stage after this age. 



ens. After the central jet is totally destroyed at a certain 
stage of FRII evolution, the hotspot vanishes at the same 
time. This is in agreement with observations which indicate 
that the hotspot luminosity decrea ses with the linear size of 
the FR II source (jPerucho fc Martiii200a ') . When the hotspot 
vanishes, the object will cease to be a "proper FRII and will 
most likely resemble a lobed FRI. A sketch of this process 
is described in Figure [T] 

Most analytical FR II models suggest the lobe is formed 
by the p articles inj ected from the central jet through the 
hotspot (|Falle|[l99ll ). Although numerical simulations show 
that there is efficient mixing between the lobe and the 
shocked ambient medium fe.g. .Scheck et al..i2002 ). their sim- 
ulations are not long enough and only represent a relatively 
early stage of the lobes. X-ray inverse Compton measure- 
ments of FRII radio lobes show that the strength of the 
magnetic field in the lobe is close to the equipartition value, 
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which suggests that the FRII lobes do not contain an ener- 
getic a lly dominant proton population l|Kataoka fc Stawar j 
I2OO5I : ICroston et all I2OO5I ). Meanwhile, the lobe internal 
pressures are in good agreement with the environment pres- 
sures, which suggests that there is no need for substan- 
tial m ixing to provide the required pressures (jBelsole et alJ 
l2007t ). These are the evidence showing that the interactions 
between the FRII lobes and their environments on large 
scales are not significant. Therefore, as our model here is 
based on previous analytical models and we are only consid- 
ering FR lis at a late stage of their evolution, we neglect the 
mixing and assume all the lobe material is from the central 
relativistic jet. As the lobe occupies a much bigger space, 
the density in the lobe is much lower than the density of 
central relativistic jet. 

The AGN active time is thought to be around a few 
10* yr, and the maximum size of FR II objects is observed 
to be a few Mpc. In order to decide whether the interaction 
between the jet and the lobe can ever be a significant factor 
in the evolution of an FRII, we therefore need to consider if 
entrainment could conceivably destroy the central jet on this 
time and/or length scale. The mixing-layer model for FRI 
objects developed by W09 discussed the interaction between 
a laminar jet and its environment. It also predicts the po- 
sition where the laminar jet disappears. We can therefore 
apply the same model to the FRII case and study the in- 
teraction between a jet and its lobe in the relativistic limit. 
This allows us to place interesting limits on the maximum 
sizes of FRII sources due to entrainment. In this section, 
we will first outline the analytical FR II lobe model and the 
W09 mixing-layer model, and then present the results with 
typical values for the parameter. 

2.1 The self-similar model for FRII lobes 



KA97 and iKaiser et all (| 19971 ) have established a success- 
ful model for FRII radio sources describing their dynamics 
and evolutions. In this section, we summarise the important 
features of the model. 

KA97 foll ow t he basic dyn amical picture proposed by 
IScheueij (|1974 ) and lFalld (|l99ll ). assuming that the laminar 
jets will end in strong shocks where the electrons are accel- 
erated. The electrons pass through the shocks and subse- 
quently inflate a lobe with a uniformly distributed pressure. 
The jet is in pressure-equilibrium with the lobe, and KA97 
showed that the lobe then expands in a self-similar way. 

The evolution of the lobe size is determined by a bal- 
ance of the ram pressure of the lobe material and that of the 
medium surrounding the host galaxy, which is pushed aside 
by the jet. The density distribution outside the core radius, 
fflo, is approximated by a power-law, p{x) — poaQX~" , where 
X is the radial distance from the central AGN and po is the 
density in the core radius, ao- KA97 suggested that, for typ- 
ical radio galaxies, po — 7.2 x 10~^^ kgm~'^ at ao — 2kpc. 
These values may vary for different sources, but we will later 
show that the precise numbers here are not important in 
our model. The exponent < a ^ 2 is constrained by both 
theories and observations. X-ray observatio ns find that the 
exponents for most cluste rs are close to 1.5 (IVikhlinin et al.l 
l2006l : ICroston et al.|[20oi '). so we adopt q = 1.5 for the mo- 
ment and will discuss the effects of adopting different values 
later in Section 3.2.1. 



Having set the density profile above, we can express the 
length of the lobe by: 



(1) 



where Qo is the jet power, A — poa^ , t is the jet age and ci 
is a constant given by equation (25) in KA97. The pressure 
of the lobe also evolves with the jet age and can be written 



Pc = 



18cf 



(F, + 1)(5- 



A-- 



(2) 



Fx is the adiabatic index of the external medium, which 
is set to be 5/3 here. Rt is the axial ratio, which is nor- 
mall y distributed between 1.3 and 6, with an average value 
of 2 l|Leahv fc Williamslll984 ). For simplicity, we adopt this 
value initially and will discuss the effect of adopting different 
values in Section 3.2.2. 

As the jet grows in a self-similar way, we can express 
the volume of the lobe hy V = ■kL^/{4R^). The particles 
injected into the jet are believed to be highly relativistic, 
so the rest mass injected into the lobe is given by mo — 
(5oi/(c^7j), where 7j is the Lorentz factor of the particles 
injected. With the expression of V and mo, the density in 
the lobe is given by: 
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V 
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t' 



(3) 



We refer the reader to KA97 for a detailed deriva- 
tion/explanation of the equations above. 



2.2 Entrainment and the mixing-layer model 

W09 model FRI sources with a mixing-layer structure in 
which a laminar jet interacts with its environment by form- 
ing a turbulent mixing layer at the interface between the 
two regions. This growing shear layer continuously entrains 
and mixes material from the jet and its environment, un- 
til finally the laminar core has been completely eroded and 
disappears. The structure of the different layers is then de- 
termined by using relativistic fiuid mechanics and applying 
the relativistic conservation laws of mass, momentum and 
energy. 

In this paper, we borrow this basic picture to estimate 
under what conditions the central jet of an FR II object may 
disappear. We assume that an FRII object evolves follow- 
ing KA97 model. Its central jet is therefore embedded inside 
the lobe and is presumably subject to entrainment from the 
lobe. The interaction of the jet with the lobe, and its sub- 
sequent evolution, are described by the relativistic mixing- 
layer model from W09. It is important to note that in the 
case of FRII type objects, the central jets are not in direct 
contact with the environment, in contrast to the model for 
FRIs presented in W09. Unlike the external medium, the 
properties of the material inside the lobe, e.g. the pressure 
and the density, are assumed to have uniform distributions. 
Thus we adopt the constant environment case of the W09 
FR I model with pe = Pc and p = pc, where pc and pc are 
given by Equations (2) and (3) respectively. 

For this simplified case, the three relativistic conserva- 
tion laws are re-written as: 
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Vs ~ 1 



(4) 



(5) 



ltPs{rt[x)~rj{x))-F{x). (6) 



Based on the equations above, the radius of the central 
jet, Tj could be expressed as a function of the distance away 
from the central AGN, x: 



2 



^(^)(r, 



(7) 



where Fj = Fs =4/3 are the adiabatic indices of the mate- 
rial inside the central jet and the shear layer, ro is the initial 
radius of the jet at the brightening point, which we assume 
is a constant equal to 100 pc throughout the life time of the 
jet. 

is defined as the ratio between rest mass energy and 
non-relativistic enthalpy for jet material. In principle, as jet 
pressure decreases with jet age, the value of should in- 
crease. However, this value may vary for different sources 
and it is hard to estimate from the observation. W09 ob- 
tained — 13.4 by applying this entrainment model to 
3C31, which is an old FRI source. Considering we are dis- 
cussing FR II sources at the late stage of their evolution, we 
assume a common value of — 10 for simplicity in our 
calculations. 

F{x) — cg{x) /[■Kp{x)] is defined in W09, where g{x) = 
Jg pvent • ndS is the mass entrainment function (n is the 
normal direction of the unit surface dS). Taking Equations 
(2) and (3), we find that F[x) is is given by: 



F{x) 
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As discussed in W09, entrainment is mainly due to turbu- 
lent motions, so the entrainment velocity Dent is closely 
related to the sound speed, Cc in the lobe. This is a con- 
stant throughout the jet lifetime, as the lobe is undergoing 
adiabatic expansion. We set Vent ~ r]Cc, and rj is the en- 
trainment efficiency which is set to be 0.5 here. KA97 and 
CR91 defined and used an entrainment efficiency in much 
the same way and argued for an upper limit of 77 < 0.26 in 
their non-relativistic mixing layer model. Our model is built 
under different conditions, and we choose a default value of 
?7 = 0.5. We discuss this issue in more detail in Section 3.2.3. 
As Jg dS is also a function of rj{x), we find that rj{x) is just 
a function of 7^ and could be solved numerically. This in- 
terestingly shows that Lmax, the maximum distance that an 
FRIT object can reach, where rj{Lmax) ~ 0, only depends 
on the Lorentz factor, 7j. 

13 = v/c and 7 = (1 — are measures of the bulk 

velocity. and jSs refer to the velocities in the central jet 
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Figure 2. The maximum length of FRII jet as a function of Pj, 
with a = 1.5 and ijy = 2. 



and the shear layer defined by W09 respectively. The analy- 
sis of some typical FR I sources indicate that bulk velocities 
are /3 ~ 0.8 - 0.9 where the jets first brighten abruptly and 
decelerate rapidly to speeds of /9 ~ 0.1 - 0.4 where recoUima- 
tion takes place fe.g. | Laing fc Bridlell2002al:ICanvin fc Laind 



I2OO4I : ICanvin et aij|2005l : iLaing et al.ll200g ). Both the ana 
lytical models and numerical simulations suggest that the 
jet velocity has a transverse structure. However, it has a 
restricted range and d oes not evolve significantly with x 
(|Laing fc Bridle! l2002bl ). so it is reasonable to use a quasi- 
one-dimensional analysis and adopt velocity values averaged 
across the jet cross-section in each region (WOO). In this pa- 
per, we adopt fis = 0.3 as a typical value. The value of j3j, 
which is a key value deciding the maximum length of FR II 
object, will be discussed in the following section. 



3 THE MAXIMUM SIZE OF FRII JETS 
3.1 Results 

The model discussed in the last section contains various pa- 
rameters and most of them have been fixed based on pre- 
vious observations or theoretical work. We first concentrate 
on the maximum length that an FRII jet can reach with 
given jet bulk velocity. Figure [2] shows the distribution of 
Lmax as a function of 7^. The diagram suggests that, if the 
particles in the central jet are slow, Lmax is small, and the 
FRII structure will be easily destroyed by entrainment at 
an early stage of its evolution. For example, an FR II with 
7j — 2 (Pj ~ 0.87) can survive only to a maximum length 
of ~ 5 kpc. As f5j increases, the jet can survive longer and 
the maximum jet length increases. An FRII with 7^ — 15 
(Pj ~ 0.998) can reach as far as 3000 kpc. This upper limit is 
sufficient to cover almost all the FRIIs currently observed. 

Relativistic beaming has been observed for many ra- 
dio sources, and it is widely accepted that powerful FRIIs 
have a very high bulk velocity. Direct measurement of jet 
bulk velocity is difficult, but several attempts have been 
made. iHough et al.l l|2002l ) observed the parsec-scale regions 
of 25 quasars (with FRII morphology) in the 3CRR sam- 
ple and estimated their bulk velocities to fall in the range 
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Figure 3. The distributions of 3CRR and 7CRS sources on the 
linear size-redshift plane. The plus signs are 3CRR sources and 
the square signs are 7CRS sources. 



of 7, » 5 — 10. X-ray obs ervations d o ne bylSambruna et al.l 
require a 7j ~ 10. iHardcastld (|2006l ) applied beamed 
inverse-Compton model to a sample of X-ray jets and shows 
the require d bulk velocity on pars ec scale could be as high 
as 7j ^ 15. iBeeelman et al.l l|2008l ) obtained similar result s 
implying 7j ^ 2 and possibly as high as ~ 50. |jorstadl (120081 ) 
studied a sample of radio loud galaxies and stated that the 
mean 7j is about 5 for radio galaxies, about 13 for BL Lacs, 
and about 20 for quasars. However, all these estimations of 
the bulk velocities are quite rough and model dependent. For 
example, although iHoueh et alT(|2002l ) measured the appar- 
ent velocities for individual sources, their orientation angles 
were not well constrained. Hardcastle (2006) used a fixed an- 
gle for all their sources. Considering that the apparent veloc- 
ities of all sources are superluminal, the value of orientation 
angle plays an important role on calculating the real bulk 
velocity. The uncertainty of the orientation angles together 
with the uncertainty of apparent velocity measurements im- 
ply huge uncertainties on the calculated real jet velocities 
and real jet sizes. We therefore cannot convincingly test our 
model by comparing predictions with observations of indi- 
vidual systems. Nevertheless, we can check for the overall 
consistency in a statistical sense here. Observational sam- 
ples of FRlIs (e.g. 3CRR, 6CE and 7CRS) show that most 
objects have sizes between lOkpc and 1000 kpc. These size 
scales refer to a 7^ > 2 — 10 based on our model. The largest 
object in the 3CRR sample has a length of ~ 4000kpc, which 
refers to > 17. All these values of Lorentz factor required 
can be fulfilled by the obs e rvatio nal results described above. 

iMullin fc Hardcastle! (|2009l ) analysed a complete sam- 
ple of FRUs and suggested that although the jet bulk 
speed on parsec scales can be as high as 7 > 10, it might 
be much smaller on kpc scales. Their fitted model gave a 
Lorentz factor around 1.18 — 1.49, corresponding to a speed 
of 0.53c — 0.74c. This is also consistent with some earlier 
work implying that the large-scale jets only have moder- 
ately relativistic bulk speeds (e.g. Wardle fc AaronI 1 19971 : 
iHardcastle et al" 1 I1999I : lArshakian fc LongaiJ |2004| ). How- 
ever, this deceleration phenomenon can be naturally ex- 
plained by our model here: the speeds observed at kpc scale 



may refer to the region containing both the central jet and 
the mixing layer. As the bulk velocity in the layer can be 0.3c 
or even lower, it is possible for us to observe a moderately 
relativistic velocity when the boundary layer dominates the 
jet on kpc scale. 

We can also check the maximum age of an FRII jet. 
With the observational constrained Lorentz factors, our 
model predicts that most FRIIs can have a maximum size 
between 10 kpc and a few 1000 kpc. lOwsianik fc Conwavl 
( 1998) observed young and powerful Compact Symmetric 
Objects, (which are believed to be the progenitors of ra- 
dio galaxies,) and obtained a hotspot advance speed of 
~ 0.03c— 0.3c at a scale of ~ 100 pc. Deceleration of the head 
may happen as a slower spe ed of '-^ 0.0 2e was obtained for 
kpc scale jet Cygnus A (Car illi fc Barth el 1996). Therefore, 
we can estimate a maximum age range from several 10^ yr 
up to 10* yr, which indicates that an FR II jet can be de- 
stroyed anytime between these time scales, but it can hardly 
survive beyond 10* yr. This range of maximum age agrees 
well with other analytical work, e.g. iBird et aL I r2008), who 
found an average jet age of 1.2 x 10^ yr. 

The most interesting conclusion from our model is that, 
as most parameters of the model are restricted in a small 
range by the observations, the maximum length of an FRII 
object is mainly determined by the bulk velocity of the cen- 
tral jet. Jets with a similar bulk velocity will have a similar 
Lmax- Here we would like to investigate two complete sam- 
ples, 3CRR and 7CRS, to see if their source distributions 
support this idea. It is hard to obtain accurate jet proper- 
ties for individuals, but as complete samples, we can assume 
the overall distributions of the jet properties (e.g. jet size, 
age and envir onment density) r eprese nt the average value at 
each redshift. IWang fc Kaiseil (|2008l ) find there is a strong 
relation between the redshift and FR II environment density 
(A oc (1 -I- z)^ *). If Lmax IS affected by the environment, we 
should observe a strong relation between the average jet size 
and the redshift. Figure [3] shows that the overall jet size only 
slightly decreases by a factor of ~ 3 from local universe to 
z — 2. This small decrease may be more due to the selec- 
tion effect: The sources are fainter when they get older and 
larger, so at high redshift, the sources easily fall below the 
flux limit before they grow to a fairly large size. The obser- 
vations are not in contrast with our model's prediction, but 
better obtained jet properties or larger complete samples are 
surely desired in the future work to provide a better support 
for our model. 

Our model also indicates that L,nax does not depend 
directly on the jet power. Although the jet power is a func- 
tion of 7j, it is also determined by the inject rate of rest 
mass. Again, we employ the 3CRR and 7CRS samples. As 
discussed in the last paragraph, we assume that at each red- 
shift, 3CRR and 7CRS sources have similar age and en- 
vironment density distributions as they are both complete 
samples. Therefore, as the 7CRS sample has a much lower 
flux limit, the sources in the 7CRS sample are generally less 
powerful than that in the 3CRR sample. In Figure O we 
can not see a significant difference between the size distri- 
butions from the two samples, which is in agreement with 
our prediction. However current observational samples are 
small with poor statistics. Larger samples with lower flux 
limit need to be considered in the future work. 
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Figure 4. The same diagram with Fig. [2] but with different value 
of a {Rj' = 2). The dotted, solid, dashed, dash-dotted line refer 
to o = 1.9, 1.5, 0.7, 0.0 respectively. 



3.2 Dependencies on parameters 

In the previous section, we discussed about how the maxi- 
mum length of a FRII object depends on the bulk velocity 
of the central jet. However, Lmax may also depend on other 
parameters which we set to be constants in Section 2. Some 
of the parameters have not been well constrained by observa- 
tions, (e.g. To and so we leave them as constants in this 
paper. Ifowever, a and Rt are well constrained and studied 
for radio jet evolution. Meanwhile the value of ry has been 
discussed a lot in previous work. Therefore, we will focus on 
these three parameters and discuss how their values could 
affect the calculated Lmax- 



3.2.1 The power-law index of environment density 
distribution, a 

iFalld i 199ll ) has shown that for a > 2, a jet shock could 
not form. X-ray observations confirm that the value of a 
should be between and 2, but the values for individual 
objects may vary significantly. In this context, we also need 
to consider the core radius, ao, inside which we take the jet 
to be surrounded by a c onstant density environment (a = 
0). ICroston et al.l (l2008l ) find that some FRI sources have 
fairly flat environments up to 100 kpc. Although this may 
not apply to FRIIs, it is still important to investigate how 
sensitive our results are to the adopted value of a. 

In order to answer this question, we calculate Lmax as a 
function of 7j for different values of a and plot the results in 
Figure |4l From this, we find that a has only a very small ef- 
fect on the relation between Lmax and 7j . If the jet is located 
in a flatter environment, it will have a slightly bigger max- 
imum length and vise versa. This result is in line with our 
assumption that the maximum length of an FR II does not 
depend directly on its environment properties. Currently, 
there are not many sources with well-determined environ- 
ment density profiles, but in the future it will be well-worth 
checking if a strong relation between a and average jet size 
exists, since this would clearly contradict with our model. 
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Figure 5. The same diagram with Fig. (2] but with different value 
of Rj- [a = 1.5). The dotted, solid, dashed, dash-dotted line refer 
to Rt = 1.3, 2.0, 4.0, 6.0 respectively. 

3.2.2 The axial ratio, Rt 

We set Rt = 2 as a constant in Section 2, but it may also 
have different values for different objects and lead to dif- 
ferent Lmax as well. Our evolutionary model is self-similar, 
so we assume Rt is constant throughout the jet lifetime. 
However, this may not be true in detail, especially during 
the late stages of FR II evoluti on. Some sim ulations sug- 
gest that Rt changes with time (|Krausell2005h . and there is 
also stron g observational evi dence that larger sources have 
larger Rt (|Mullin et al.|[200^ ). It is therefore also important 
to check how sensitive Lmax is to a variable Rt. 

Figure [5] shows the Lmax — 7j diagram with different 
values of Rt. It shows that a fatter jet should have a larger 
maximum length. As Rt is used for calculating the prop- 
erties inside the lobe, which directly associate with entrain- 
ment process, it is reasonable to have a bigger influence on 
Lmax than a. However, the minimum and maximum sug- 
gested values of Rt only change the final Lmax by a factor 
of around 3, so our assumption of constant Rt should still 
be a tolerable approximation. 

3.2.3 The entramment efficiency, rj 

The entrainment efficiency is another important parameter 
constraining the maximum size of an FRII jet, but it is 
difficult to obtain from either theoretical analysis or obser- 
vations. CR91 discussed the entrainment efficiency in the 
context of their model and showed that this is limited either 
by the ability of the environment /jet to supply material for 
the mixing layer, or by the maximum possible growth rate 
of mixing layer itself. They defined three different regimes, 
each associated with a particular upper limit on the entrain- 
ment efficiency. In our model, the environment is actually 
the cocoon and pc <S Pj, so the appropriate regime is the 
environment-limited regime. As the system is in a pressure 
equilibrium, we must then have a sound speed in the cocoon 
(environment) much larger than that in the jet. Therefore, 
we can only set a common upper limit of < 1. The actual 
value of rj for individual sources is more difficult to gen- 
eralise as a number of instabilities (e.g. Kelvin-Helmholtz, 
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current-driven and so on) together with jet properties can 
play important roles in setting the value of rj. As a result, 
in our actual calculations in this paper, we have simply as- 
sumed a typical value of = 0.5. Different values of in 
the allowed range will result in a large Lmax coverage, so we 
do not intend to plot a Lmax — 7j diagram here. However, 
it is easy to see that a lower (higher) entrainment efficiency 
will lead to a slower (faster) destruction of the central jet 
and to a larger (smaller) maximum jet size. For example, if 
r] = 0.9, a jet need to have 7j > 12 to reach a maximum size 
of 1000 kpc, but for rj — 0.1, the same Lmax can be reached 
with only > 5. 

3.3 Comparison with previous work 

KA97 also discussed the stability of jet based on the mixing 
layer model from CR91. Given that the W09 model used 
here is also based on CR91, it is interesting to compare our 
results to those obtained by KA97. 

KA97 found that in an environment with a — 0, the jet 
can easily be destroyed and can only be stable up to 2.6 kpc. 
With 7j = 2 (KA97's default assumption), our model pre- 
dicts a slightly lager Lmax of around 5 kpc for this a. Ifow- 
ever, KA97 also argued that the Lmax will increase with a, 
i.e. jets in the environment with a steeper density gradient 
can grow to larger sizes before becoming unstable. This con- 
trasts with our model, in which the value of a plays only a 
minor role in setting Lmax, with being the key parameter 
instead. 

There are two main differences between the models 
from KA97 and this paper. First, KA97 adopt mixing-layer- 
limited regime of CR91 model, and obtain an upper limit 
of r; < 0.26 for the entrainment efficiency. However, as 
we discussed in the last section, the environment-limited 
regime should be more appropriate here. Meanwhile, the 
original CR91 model applied by KA97 was designed for non- 
relativistic cases, whereas the W09 model is based on the rel- 
ativistic conservation laws. In the non-relativistic case, the 
energy and momentum are mainly decided by the density 
of the lobe, which strongly depends on lobe volume and a. 
However, in the relativistic case, the momentum and energy 
are dominated by the relativistic component, so the Lorentz 
factor is the key parameter. 

A number of numerical simulations have studied jet in- 
stabilities in more detail taking into account more factors 
other tha n the Lorentz factor a s considered in our model. For 
example, iMizuno et all (|2007l ) suggested that the distribu- 
tion of magnetic field is crucial for determining the K elvin- 
Helmholtz stabilisa tion. Moreover, iRossi et al.l (|2008l ) and 
iMeliani et al.l (|2008l ) found that the environment/jet density 
contrast is important in determining the instability evolu- 
tion and entrainment properties. However, all these simula- 
tions only represent the early stage of jet evolution. At this 
stage, the lobe is still small and not well established. Mean- 
while, the mixing between the environment and the lobe may 
be significant, and the lobe density is much higher than that 
at the later stages of jet evolution, which are considered in 
this paper here. Therefore t he properties of the e nvironment 
play more important roles. IPerucho et al.l (|2004l ) performed 
long-term simulations for the relativistic jets, and consid- 
ered a slightly overdensed environment and lobe, which is 
more closed to our case here. They also found that the jet 



Lorentz factor is an important parameter deciding the non- 
linear stabilities of the jets. With the same thermodynami- 
cal properties, the jets with smaller Lorentz factors start to 
mix and transfer momentum to the environment at a earlier 
stage. However, as they tried various models with different 
values of thermodynamical properties, they claimed that a 
number of other parameters (e.g., jet-to-ambient enthalpy 
ratio, temperature and jet internal energy,) also affect the 
jet stabilities and long-term morphologies. Although in our 
model here, we only consider the jet Lorentz factor as the 
key parameter, it is worth reiterating one particular simpli- 
fication we make, which is that we attribute all kinds of jet 
instability evolutions to the growing of the shear layer be- 
tween the central jet and the lobe. In reality, the evolution 
and the growing speed of the shear layer must surely be a 
function of several other physical para meters an d may there- 
fore vary in time and between sources (Peruclio et al. 2005i). 
We make this simplification here purely to keep the model 
analytically tractable, although with sufficient data it may 
become possible to reconstruct the dependence of instabil- 
ity evolutions on various jet parameters from observations 
in the future. 



4 THE EVOLUTION FROM FR II INTO FR I 
SOURCES 

The existence of a maximum size for FRII sources due to 
the erosion of their laminar jets raises an obvious ques- 
tion: what happens to an FRII object that reaches this 
limit? iBicknelll (f 19941 ) suggested that the transition between 
FRI/II is due to the transition from subrelativistic to rela- 
tivistic flow caused by entrainment. As the death of FRIIs 
in our model here is also due to the entrainment, we argue 
that the FRIIs reaching their maximum sizes are likely to 
evolve into FRIs. In this section, we will outline a simple, 
but plausible scenario for the transition of a radio galaxy 
with an FRII morphology to one with an FRI morphology. 
The basic idea is sketched in Figure [T] 

When a stable radio outflow is born at time to, it ex- 
hibits an FRII structure with a laminar flow embedded in- 
side a lobe and a hotspot at the end. At stage ti, when the 
jet length Lj is smaller than the maximum length Lmax, the 
outflow grows with age, following the KA97 picture. Mean- 
while, however, the central jet continuously suffers entrain- 
ment from the lobe, and the structure of the outflow can 
be described by our model here. The outflow evolves with 
an FRII morphology until it reaches Lmax at age tmax- At 
this time, the central jet is totally eroded, and the hotspot 
vanishes. The detailed evolution process of the radio outflow 
at this stage is described in Section 2. 

When the outflow evolves to an age of t2, where t2 > 
tmax, the shear layer dominates the end region of the jet and 
the hotspot vanishes. A weaker shock and the lobe structure 
may still exist with plasma injected into the lobe after the 
shock from the end of the jet. The expected structure at this 
stage is reminiscent of a typical lobed FRI source. As the 
outflow becomes even older, the energy from the shear layer 
can hardly support the lobe structure or the working surface 
of the shock at the end of the jet, so the plasma will form a 
turbulent tail with the lobe disappearing either because it is 
reflUed from the environment or because it simply runs out 
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of energy. At the end of this evolution stage, we will observe 
a naked tailed jet like 3C 31. The jet is in direct contact with 
the environment, and a mixing shear layer is formed. At the 
same time, the laminar part may shrink again as the density 
of the environment is higher than that of the lobe. 

Please note that we are not claiming that our transi- 
tion scenario here is the only way to generate FRIs. The 
precursors of FRIs ma y also include weak CSS sources 
Perucho fc Martil l2007l l. ets hitting dense environments 



Meliani et al.l 2008h , and weak FR lis reaching the pressure 



equilibrium with the ir environment inside the core region 
ijKaiser fc Best|[2007^ . These works are not in contrast with 
our work here as we are considering if FR lis can evolve into 
FR Is at the late stage of their evolution. Our transition sce- 
nario is just complementary to all the work above, providing 
a new plausible way for powerful FRIIs to develop into FRIs 
later in their lives. More studies are still needed in order to 
fully understand the FRI/II dichotomy. 



5 CONCLUSION 

We have embedded a mixing-layer model originally devel- 
oped for modelling FRI jets into a self-similar model for 
FR II radio lobes to study the effect of entrainment on the 
central jets in FR II objects. We find that, for reasonable 
parameters, the growing mixing layer between the central 
jet and the radio lobe could play an important role during 
the evolution of FR II objects. 

The maximum length that a jet can reach is decided 
mainly by the bulk velocity of the particles in the central 
jet, Pj, but not directly related to the environment or the 
jet power. We find a maximum length of ~1000 kpc for 
Pj — 0.997, assuming an entrainment efficiency of 77 = 0.5, 
an environment index of a = 1.5 and a lobe axial ratio of 
Rt = 2. If the jet is located in a flatter environment with a 
smaller a or the jet is fatter with a smaller Rt, its maximum 
length will be larger. 

We have also sketched the likely evolution of FRII 
sources after they reach their maximum size. Once the 
hotspots are extinguished, such sources will initially look like 
lobed FRI objects. However, ultimately their lobes will stop 
getting fed from the jet, become turbulent and be refilled by 
the environment. At this point they will emerge as classic, 
SCSI-like FRI sources. This simple scenario suggests a new 
evolutionary connection between FR I and FR II sources and 
may help to shed new light on the FRI/II dichotomy. 

The FRI/II transition process suggested by our model 
may not depend on the environment directly, but there is 
some observational evidence showing t hat FRIs and FRIIs 
may inhabit different in environments (jPrestaee fc Peacockl 
Il988l ). Thus we cannot totally rule out the influence from 
the environment. There might be relations between envi- 
ronment properties and model dependent parameters (e.g. 
7j), so that the environment can affect Lmax indirectly. This 
is an interesting open question for the future work once we 
have more observational data for the jet environment. 

In closing, we stress that the picture we have developed 
here - especially that of the evolution beyond tmax - is still 
basically a toy model. In the future work, we are planning 
to model the evolution of the jet from FR II to lobed FR I to 
tailed FR I in more detail. Our final goal is to build a unified 



model for all types of radio galaxies and track how they 
evolve and morph into each other across the P-D diagram. 
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